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This paper is concerned with the long wavelength instabilities (infrared catastrophes) occurring in 
Bose-Einstein condensates (BECs). We examine the modulational instability in "cigar-shaped" (ID) 
attractive BECs and the transverse instability of dark solitons in "pancake" (2D) repulsive BECs. 
We suggest mechanisms, and give explicit estimates, on how to "engineer" the trapping conditions 
of the condensate to avoid such instabilities: the main result being that a tight enough trapping 
potential suppresses the instabilities present in the homogeneous limit. We compare the obtained 
estimates with numerical results and we highlight the relevant regimes of dynamical behavior. 



I. INTRODUCTION AND SETUP 

Infrared catastrophes, or long- wavelength instabilities (LWI) as they are otherwise known, are ubiquitous in physical 
phenomena. From Magneto-Hydro-Dynamics [1] to chemical models [2], and from quantum systems [3] to fluid 
mechanics [4], polymer physics [5] and plasmas [6], large scale modulations may destabilize the system of interest. 
On the other hand, many of these instabilities (and their thresholds) have been quantified in the context of nonlinear 
models. 

In the past few years, another context that can be accurately modeled by nonlinear partial differential equations 
that are known to possess LWI, has become experimentally tractable. This setting is, in particular, the one of Bose- 
Einstein condensates (BECs) [7] whose experimental realization has led to an explosion of interest in the field of 
atomic matter waves and of nonlinear excitations in them. Such nonlinear waves have been recently experimentally 
generated in BECs, namely dark [8] and bright [9,10] solitons. Also, two dimensional excitations, such as vortices [11] 
and lattice patterns thereof [12] have also been obtained experimentally, while other nonlinear waves, such as Faraday 
waves [13], ring solitons and vortex necklaces [14] have been theoretically predicted. It is worth noting here that the 
nature of Bose-Einstein nonlinear matter waves depends crucially on the type of the interatomic interactions: dark 
(bright) solitons can be created in BECs with repulsive (attractive) interatomic interactions, resulting from positive 
(negative) scattering length. 

In the present work, our scope is to re-examine the LWI in the context of Bose-Einstein condensates and exploit 
the consequences on the LWI of the inhomogeneity induced by the trapping potential. It is well-known [7] in this 
setting that the condensates are formed under appropriate confining conditions that typically consist of magnetic 
traps modeled by parabolic potentials. Our scope is then two-fold. On the one hand, we aim at illustrating the 
potential for the instabilities and at examining their dynamical development when they exist. On the other hand, 
the magnetic confinement provides an additional "trapping" length scale to the problem whose competition with the 
instability length scale may disallow the dynamical manifestation of the LWI. In this way, we will propose how to 
engineer trapping conditions so as to prevent unstable dynamical evolution. 

We will examine two benchmark examples of long-wavelength instabilities in BECs: the modulational instability 
occurring in attractive, one-dimensional (ID) BECs and giving rise to the formation of bright matter-wave solitons, 
and the transverse ("snaking") instability of dark soliton stripes formed in repulsive two-dimensional (2D) BECs. It 
should be noted here that the genuinely three-dimensional (3D) condensate can be considered as approximately ID 
if the nonlinear inter-atomic interaction is weak relative to the trapping potential force in the transverse directions; 
then, the transverse size of the condensates is much smaller than their length, i.e., the BEC is "cigar-shaped" and can 
be effectively described by ID models [15,16]. Similarly, if the transverse confinement is strong along one direction 
and weak along the others, then for this "pancake-shaped" BEC, 2D model equations are relevant [17]. 

Close to zero temperature, it is well-known that the 3D Gross-Pitaevskii (GP) equation [7] accurately captures the 
dynamics of the condensate. For cigar-shaped BECs, the model equation is effectively ID and can be expressed in 
the following dimcnsionless form, 
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l - = --— + a\ufu + V( X )u, (1) 

where u is the macroscopic wave function (normalized to 1), and a = a/\a\ = ±1 is the renormalized scattering length, 
which is positive (negative) for repulsive (attractive) condensates. In this equation, t and x are respectively measured 



1 



in units of 9/16e 2 u>± and 3a^/4e, where ujj_ is the confining frequency in the transverse direction, aj_ = yjh/moj^ 
is the transverse harmonic-oscillator length and e = N\a\/a± is a small dimensionless parameter, N being the total 
number of atoms [16]. Finally, 

V(x) = ^-x 2 (2) 

where O = (9/16e 2 ) • (ui x /u)±) is the frequency of the magnetic trap potential V in our dimensionless units, and u x is 
the axial confining frequency. Similarly the 2D model for the pancake-shaped condensate assumes the form: 

.9u 1 . l9 Ir , . 

i— = --Au + a\ufu + V(r)u (3) 

where now the role of the axial frequency /spatial variable is played by the radial one (r = y/ 'x 2 + y 2 ) and the 
normalized variables are connected to the dimensional ones similarly to the ID case, but with the role of lo x (ui±) now 
played by uo±_ (io z ) . 

We structure our presentation as follows: in Section II we study the suppression of the modulational instability 
for ID cigar-shaped BECs in presence of a tight confining magnetic potential, while in Section III we investigate the 
transverse instability for 2D trapped BECs. In both cases we provide a criterion (in terms of the trap parameters) 
giving the conditions under which LWI can be avoided. Section IV is devoted to our concluding remarks. 



II. MODULATIONAL INSTABILITY 



It is well-known (see, e.g., the recent work [18] and references therein) that, in the absence of external potential, 
the continuous-wave (cw) solution u = w cxp(— icm 2 ^) of amplitude uq, of Eq. (1) becomes modulationally unstable 
when perturbations of wavenumber k < k cr = 2|a| 1 / 2 u are imposed. This can be equivalently interpreted as follows: 
when length scales 

A > X cr = — ^= (4) 



become "available" to the system, then the modulation over these scales leads the solution to instability. 

However, in the presence of the magnetic trap, there is a characteristic scale set by the trap, namely the BEC axial 
size, XbeCi which depends on the trapping frequency Q. When Xbec < X cr suppression of the modulational instability 
is expected. To estimate Xbec i n & specific setup, we will examine a protocol relevant to the recent experiments 
conducted by the Rice [9] and Paris [10] groups, which has stimulated a considerable amount of theoretical attention 
[19]. In particular, we start with a ID repulsive condensate, whose ground state wavefunction is approximately in 
the so-called Thomas- Fermi (TF) regime [7], and subsequently change the interaction into an attractive one. This 
is experimentally achieved using the so-called Feshbach resonance: an external magnetic field is used to modify the 
scattering length of the interatomic interactions [20]. In our example, we use as initial condition the ground state 
of Eq. (1) with a = 1, which is in the TF approximation \u\ 2 w fi — V(x); then at t — we change the sign of the 
scattering length, i.e., we set a = —1. Therefore in this situation Xbec ~ 2y/2ji/Q. From the normalization condition 
J dx\u\ 2 = 1 one gets 

/ \ 2 / 3 

"(3) • 

The condition for the suppression of the modulational instability Xbec < X cr [where X cr is given by Eq. (4)] gives 
fl > 2 3 / 2 (\a\(i) 1 / 2 u /ir. In this context, the amplitude u in Eq. (6) can be well approximated as u ~ y/Ji, since the 
TF approximation is most accurate close to the center of the condensate. Therefore if ft > Q cr , where 

Her = ^3 « 0.2, (6) 

then the trapping conditions are "engineered" in such a way that the modulational instability cannot manifest itself: 
for a tight enough trapping potential the modulational instability does not occur. This is one of the key aims of this 
paper, namely to quantitatively highlight how the infrared catastrophes can be avoided in the presence of sufficiently 
"tight" trapping of the condensate. We remark that if one evaluates Xbec by using for the ground state of Eq. (1) 
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(with a = 1) a gaussian with variationally determined width, one has an estimate of f2 cr in good agreement with Eq. 
(6). We also notice that the critical length X cr is a few healing lengths £: indeed, when the the density grows from 
to Uq within a distance £, the quantum pressure and interaction energy terms are equal when l/(2£ 2 ) w Uq, which 
means that £ w (\/2uo) _1 and therefore X cr w 7r\/2£ = 4.44£. 

Typical experimental values for a 7 Li BEC are ^ = 2n x 5Hz, = 2ir x 500Hz, a = — 3a and A w 10 3 ; these 
yield e ~ 0.1 and £1 ~ 0.5, which should be a sufficiently large value to avoid the emergence of the modulational 
instability. 

We tested this prediction by means of direct numerical simulations shown in Fig. 1. The panels (a), (b), (c) show 
respectively the case of f2 = 0.3, 0.1, and 0.02. It is seen that the supercritical, "tight" trap f2 = 0.3 docs not allow 
the development of the instability for the condensate [Fig. 1 (a)] . The only consequence of the change of the sign 
of the scattering length is the excitation of an internal mode of the condensate, due to the fact that the TF cloud 
is not the ground-state for a = — 1. This results in nearly periodic oscillations of the width of the wavefunction 
W = [J x 2 \u\ 2 dx — (j x\u\ 2 )Y^ 2 , which we monitor as a diagnostic in our numerical examples. For il = 0.1 [Fig. 1 
(b)], the width oscillations are no longer periodic: the result of the modulational instability can be clearly discerned 
in the addition of an extra frequency to the motion of the condensate for subcritical trapping. We have identified 
this to be a rather smooth transition, becoming more pronounced close to the theoretically predicted critical point 
of fi « 0.2. We have also examined the development of the instability for weaker trapping conditions [f2 = 0.02, see 
Figs. 1 (c), (d)]. Then an additional frequency emerges, as the condensate is now cleaved in two pieces during the 
time evolution [see Fig. 1 (d)]. For even weaker trappings, we found that the larger the ratio of \BEc/X cri the more 
frequencies appear in the condensate density evolution, rendering it increasingly chaotic. Hence, we have demonstrated 
that a cascade of frequencies arises in the BEC dynamics due to the modulational instability and subsequent higher 
resonances that allow the breakup of the condensate and make its motion less regular. This mechanism should clearly 
be experimentally detectable; in fact we surmise that the observation of a single matter- wave bright soliton [10], rather 
than a soliton train [9], is due to the fact that Xbec was smaller in the former case, as the number of atoms was an 
order of magnitude fewer in the Paris experiment, while the ratio of the trapping frequencies was approximately the 
same. 



III. TRANSVERSE INSTABILITY 



One of the infrared catastrophes that occur in two spatial dimensions is the transverse instability of dark-soliton 
stripes in repulsive (a > 0) BECs. As a result, a dark-soliton undergoes a transverse snake deformation [21], causing 
the nodal plane to decay into vortex pairs. This instability has been examined in the context of BECs in [22], and 
the Bogoliubov spectrum of the dark soliton has been obtained; in this context, the relevant imaginary modes were 
identified to transfer the energy of the condensate to collective excitations parallel to the nodal plane destroying 
the configuration. However, as is highlighted in [23], "the explicit connection between the existence of imaginary 
excitations and a dynamical snake instability remains unclear" . Our scope here is to illustrate the criterion for the 
transverse instability and to test it against direct numerical simulations, exposing the possible dynamical scenarios 
and quantifying their dependence on the trapping parameters. 

In the absence of the potential, the transverse instability occurs for perturbation wavenumbers 



sin 2 (j> + u 2 sin 4> + u 4 — (2 sin <j> + u 2 ) 



1/2 

(7) 



where sin<fi is the dark-soliton velocity [21] and u is the amplitude of the homogeneous background, connected with 
the chemical potential through w 2 = /x similarly to the previous (ID) setting. In the case of stationary (black) solitons, 
of interest here, sin = 0, hence k cr = u^ 1 . On the other hand, for V(r) = £l 2 r 2 /2, a similar calculation as for the 
ID problem yields the characteristic length scale of the BEC (i.e., the diameter of the TF cloud) as Xbec ~ 2y/2jl/il. 
Then, the criterion for the suppression of the transverse instability is that the scale of the BEC is shorter than the 
minimal one for the instability. The corresponding condition reads 



n > (8) 

TTUo 

To obtain the minimum value of Jl we need to know how /i is connected with uq. As a first guess, in the absence of the 
dark soliton, one can assume m§~|U (close to the center of the BEC), which yields £1 > \/2/tt = 0.45. Hence, stronger 
trapping should "drown" the transverse instability and preserve dark soliton stripes on top of the Thomas- Fermi cloud 
(i.e., stable "dipole" solutions). Note that in terms of real physical units, the above mentioned critical value of ft may 
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correspond, e.g., to a weakly interacting 87 Rb pancake condensate, containing w 10 3 atoms, confined in a trap with 
uj r = 2tt x 5Hz and io z = 2ir x 50Hz. 

We have numerically tested this condition, finding it to be an overestimate of the critical trapping frequency for the 
transverse instability, which is £l cr w 0.31. This result shows that the numerically found Fig. 2 shows a dynamical 
evolution example of the dipolar solution, initialized with a tanh(y) imposed on the TF cloud, for £1 = 0.35 (left) 
and £1 = 0.15 (right). Both snapshots show the contour plot of the square modulus of the wave function at t = 1000 
(left) and t = 200 (right). Clearly, in the former case the transverse instability is suppressed, while in the latter 
a vortex pair has been formed demonstrating the dynamical instability of the configuration. We have monitored 
the asymptotic, long time evolution of the dipole and have observed the following interesting phenomenology: for 
0.18 < fl < 0.31, while the stripe is dynamically unstable, there is not sufficient space for the instability-induced 
vortices to fully develop; as a result, after their formation, they subsequently recombine and disappear. This behavior 
is shown in Fig. 3 (for £1 = 0.2): the two vortices formed are shown at t = 190 (left panel), while at t = 210 they 
recombine to form a transient dark stripe (right panel). This configuration is unstable and it subsequently breaks up 
to a new vortex pair (not shown here), which eventually recombines at longer times (t w 400). It is interesting to 
note that as, in the present case, the available size of the condensate is of the order of a few healing lengths £ (in our 
units, £ = 1/a/2), the two vortices formed hardly fit the condensate size (recall that the vortex core is of order of £ 
[24]). This is a possible qualitative explanation of this recombination, whose origin is the competition of the length 
scales available in the BEC. Finally, it should be noted that for £1 < 0.18, the vortices will survive in the asymptotic 
evolution of the system, and naturally the weaker the trapping the larger the number of "engulfed" vortices generated 
due to the stripe breakup. 

A question that naturally arises in the results above concerns the disparity between the critical point theoretical 
estimate for the transverse instability and the corresponding numerical finding. We believe that the disparity is 
justified by the fact that the theoretical stability analysis of [21] is performed for the infinite homogeneous medium 
(the dark-soliton pedestal), in the absence of a magnetic trap. In the presence of the trap on the one hand, the 
background is inhomogeneous, while, on the other hand, for tight traps resulting in small condensate sizes, the 
presence of the dark soliton at the BEC center, significantly modifies the maximum density. Thus, one should expect 
that the relation Uq~^l should be modified as Suq — /i, where the "rescaling" factor 8 < 1. Based on the results of 
the numerical simulations (see Figs. 2-3) we have observed that embedding the dark soliton in the BEC center reduces 
the maximum density of the TF cloud to the half of its initial value. This suggests that 8 — 1/2, which, according to 
Eq. (8), leads to a new value for £1 minimum, i.e., £1 = 1/n w 0.318. This modified criterion for the suppression of 
the transverse instability, namely £1 > 0.318 is in very good agreement with the above mentioned numerically found 
condition £1 > 0.31. 



IV. CONCLUSIONS 

We have examined case examples of long wavelength instabilities in Bose-Einstein condensates. We have revisited 
the modulational instability in quasi-lD, cigar-shaped BECs, and the transverse instability in quasi-2D, pancake- 
shaped BECs. We have advocated that trapping conditions can be engineered to avoid or induce such instabilities at 
will. Using the length scale competition of the entrainment due to trapping and of critical instability wavelength, we 
have given explicit estimates for critical values of the trapping frequency beyond which the instabilities will be absent. 
We have tested these criteria in both cases and have found good agreement with the numerical results in the ID 
case. In the 2D setting, we have explained the overestimation of the critical point, on the basis of the homogeneous 
background assumed in the theoretical estimate, as well as the modification of the value of the maximum density 
of the condensate due to the presence of the dark soliton. Our results demonstrate how to engineer the trapping 
conditions, in order to achieve supercritical regimes, devoid of long wavelength instabilities. On the other hand, for 
the subcritical regimes, we have illustrated the relevant phenomenology through direct numerical simulations and the 
cascade that leads to the instability and eventual destruction of the original coherent structure. Such results can be 
particularly useful in quantifying the selection of external conditions so as to achieve or avoid a given experimental 
outcome. 
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FIG. 1. Condensate width as a function of time for Q = 0.3 (a), f2 = 0.1 (b), and fi = 0.02 (c). In (d) we plot the atomic 
density with £1 — 0.02 at three different times: t = (solid line), t = 80 (dotted line, corresponding to point 1 in (c)) and 
t = 160 (dashed line, corresponding to point 2 in (c)). 
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FIG. 2. The panels show the contour plots of the density \u\ 2 for Q = 0.35 at t = 1000 (left) and Q = 0.15 at t = 200 (right). 
In the first case the transverse instability is clearly suppressed, while in the second it sets in, giving rise to a formation of a 
vortex pair. 




FIG. 3. Snapshots of a vortex-pair evolution in a case where snaking instability has set in (fl = 0.2). In the left panel 
(t = 190) the formed vortex pair is shown, while the right panel (t = 210) shows the recombination of the two vortices, 
resulting in the re-generation of a dark stripe structure. The latter is unstable and decays at longer times (t ~ 400). 
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